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The natural excitations of an interacting one-dimensional system at low energy are hydrodynamic 
modes of Luttinger liquid, protected by the Lorentz invariance of the linear dispersion. We show 
that beyond low energies, where quadratic dispersion reduces the symmetry to Galilean, the main 
character of the many-body excitations changes into a hierarchy: calculations of dynamic correlation 
functions for fermions (without spin) show that the spectral weights of the excitations are propor¬ 
tional to powers of 7 Z 2 /L 2 , where 72 is a length-scale related to interactions and L is the system 
length. Thus only small numbers of excitations carry the principal spectral power in representative 
regions on the energy-momentum planes. We have analysed the spectral function in detail and have 
shown that the first-level (strongest) excitations form a mode with parabolic dispersion, like that of 
a renormalised single particle. The second-level excitations produce a singular power-law line shape 
to the first-level mode and multiple power-laws at the spectral edge. We have illustrated crossover to 
Luttinger liquid at low energy by calculating the local density of state through all energy scales: from 
linear to non-linear, and to above the chemical potential energies. In order to test this model, we 
have carried out experiments to measure momentum-resolved tunnelling of electrons (fermions with 
spin) from/to a wire formed within a GaAs heterostructure. We observe well-resolved spin-charge 
separation at low energy with appreciable interaction strength and only a parabolic dispersion of 
the first-level mode at higher energies. We find structure resembling the second-level excitations, 
which dies away rapidly at high momentum in line with the theoretical predictions here. 

PACS numbers: 71.10.Pm, 03.75.Kk, 73.63.Nm, 73.90.+f 


I. INTRODUCTION 


in a manner consistent with a power law. 


Predicting the behaviour of interacting electrons is a 
significant open problem. Most progress to date has 
been made at low energies where linearisation of the 
single-particle dispersion led to construction of Fermi 
and, in one dimension, to Luttinger-liquid theories^ in 
which the natural excitations are fermionic quasiparticles 
and hydrodynamic modes respectively. The only sig¬ 
nificant progress beyond the linear approximation has 
been achieved via the heavy impurity model, for Fermi — 
and Luttinger— liquids, showcasing threshold singularit¬ 
ies drastically different from the low energy behaviour. 
In this paper we investigate one-dimensional (ID) fermi¬ 
ons beyond the linear approximation where the natural 
many-body excitations form a hierarchical structure^ in 
sharp contrast with the Fermi quasiparticles or hydro- 
dynamic modes. We obtain the dynamical structure 
factor, in addition to the already known spectral func¬ 
tion, and construct an inductive proof for calculating 
the form factors that are necessary for the dynamical 
response functions of the spinless fermion model. Exper¬ 
imentally, we demonstrate control over the interaction 
energy a ID wire manifested as a change of the ratio of 
the charge and spin velocities at low energy scales. We 
find a new structure resembling the second-level excita¬ 
tions, which dies rapidly away from the first-level mode 



Figure 1. Regions of the energy-momentum plane dominated 
by two different principal regimes of the system (bottom): 
hydrodynamic modes of the Luttinger liquid (top-right) at 
low energies (marked with cyan colour in the bottom panel) 
and the hierarchy of modes (top-left) in the rest of the plane. 
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We analyse theoretically the dynamic response func- 
tions-that probe the many-body excitations-for spin¬ 
less fermions with short-range interactions. Our ap¬ 
proach is exact diagonalisation via Bethe ansatz meth¬ 
ods: the eigenenergies are evaluated in the coordinate 
representation and form factors-for the corresponding 
eigenstates-are derived in the algebraic representation, 
via Slavnov’s formula^ On the microscopic level the ex¬ 
citations arrange themselves into a hierarchy via their 
spectral weights-given by the form factors-with different 
powers of 7 Z 2 /L 2 , where TZ is the particle-particle inter¬ 
action radius and L is the length of the system. As a 
result only small numbers of states-out of an exponen¬ 
tially large Fock space of the many-body system-carry 
the principal spectral power in representative regions on 
the energy-momentum plane, see Fig.[TJ allowing an ana¬ 
lytical evaluation of the observables. 

At small energy this hierarchy crosses over to hydro- 
dynamic behaviour, see Fig. QJ that we illustrate by cal¬ 
culating the local density of states. At low energy it 
is suppressed in a power-law fashion according to the 
Tomonaga-Luttinger theory. Away from the Fermi point 
where the Lorentz invariance is reduced to Galilean by 
the parabolicity of the spectrum, the local density of 
states is dominated by the first(leading)-level excitations 
of the hierarchy. This produces a 1/y/e van Hove sin¬ 
gularity, where 5 is the energy measured from the bot¬ 
tom of the conduction band. At even higher energies 
the second-level excitations produce another 1 / y/\s\ van 
Hove singularity on the other side of the band edge, in 
the forbidden for the non-interacting system region. 

Using this framework, we study response of the cor¬ 
related system to adding/removing a particle in detail, 
given by the spectral function. The first-level excitations 
form a parabolic dispersion, like a single particle, with a 
mass renormalised by the Luttinger parameter The 
continuous spectrum of the second-level excitations pro¬ 
duces a power-law line-shape around the first-level mode 
with a singular exponent —1. Around the spectral edges 
the second-level excitations give a power-law behaviour 
of the spectral function. For the hole edge the expo¬ 
nent calculated microscopically reproduces the predic¬ 
tion of the phenomenological heavy impurity model in 
one-dimension^ However, around the particle edge the 
second-level excitations give a power-law of a new type. 

Experimentally, momentum-resolved tunneling of elec¬ 
trons confined to a ID geometry has been used to probe 
spin-charge separation in a Luttinger liqmd i n i 12 i 38 i 39 
This separation was observed to persist far beyond 
the energy range for which the Luttinger approxima¬ 
tion is validji^ showing the need for more sophisticated 
theories^ 7 Particle-hole asymmetry has also been detec¬ 
ted in relaxation processes^ In this paper we measure 
momentum-resolved tunneling of electrons in the upper 
layer of a GaAs-AlGaAs double-quantum-well structure 
from/to a 2D electron gas in the lower layer. This set-up 
probes the spectral function for spinful fermions. We ob¬ 
serve well-resolved spin-charge separation at low energy 


with appreciable interaction strength—a distinct effect of 
the spinful generalisation of Luttinger liquid^ The ratio 
of charge and spin velocities is v c /v s ~ 1.8^ At high 
energy, in addition to the spin and charge curves, we can 
also resolve structure just above fcp that appears to be 
the edge of the second-level excitations. However, the 
amplitude decays rapidly and for higher k we find no 
sign of the higher-level excitations, implying that their 
amplitude must have become at least three orders of 
magnitude weaker than for the parabola formed by the 
first-level excitations. The picture emerging out of these 
experimental results can only be explained—though only 
qualitatively—by the hierarchy that we study for spinless 
fermions. 

The rest of the paper is organised as follows. In Sec¬ 
tion II we describe the one-dimensional model of inter¬ 
acting spinless fermions introducing a short range cut-off 
via lattice. Section III contains a procedure of finding 
the many-body eigenenergy by means of the coordinate 
Bethe ansatz. In Section IV we evaluate the form factors 
needed for the dynamical response functions. We give 
a construction of the algebraic representation of Bethe 
ansatz (Subsection IVa) and evaluate the scalar product 
in this representation (Subsection IVb). We present a 
calculation of the form factors for the spectral function 
and the dynamical structure factor for a finite chain (Sec¬ 
tion IVc). We take the limit of long wavelengths deriv¬ 
ing polynomial formulae for the form factors (Subsection 
IVd). Then, we analyse the obtained form factors estab¬ 
lishing hierarchy of excitations (Subsection IVe). Finally 
we calculate the spectral function around the spectral 
edges (Subsection IVf). In Section V we illustrate the 
crossover to Luttinger liquid at low energy by evaluating 
the local density of states at all energy scales. Section VI 
describes experiments on momentum-conserved tunnel¬ 
ling of electrons in semiconductor wires. Section VII is 
dedicated to low energies and in Section VIII we analyse 
the measurements at high energies connecting the exper¬ 
iment with theory on spinless fermions developed in this 
paper. Figures below are marked with spinless and spin¬ 
ful logos (such as those in Fig. [2]and[7J respectively) to 
indicate the structure of the paper visually. Appendix A 
contains details of the derivation of the Bethe equations 
in the algebraic representation. In Appendix B we derive 
the expectation value of the local density operator. 

II. MODEL OF SPINLESS FERMIONS 

We study theoretically the model of interacting Fermi 
particles without spin in ID, 

dx ( — — (x) Aip (x) + ULP (x) 2 J , (1) 

\ 2 m J 

where the field operators ip (x) satisfy the Fermi com¬ 
mutation relations, {ip (x ), ip^ (#')} = S (x — x'), p (x) = 
ip^ (x) ip (x) is the particle density operator, and m is 
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the bare mass of a single particle. Below we consider 
periodic boundary conditions, ip {pc + L) — ijj (x) , restrict 
ourselves to repulsive interaction U > 0 only, and take 
ft = 1. 

Non-zero matrix elements of the interaction term in 
Eq. dU) require a finite range of the potential profile for 
Fermi particles. Here, we will introduce a lattice with 
next-neighbour interaction which lattice parameter and 
interaction radius is 7 Z. The model in Eq. m becomes 


H = 


2^ (V’jV’i+i + y-V’i-i) + Upjp j+ i 


( 2 ) 


where j is the site index on the lattice, the dimensionless 
length of the system is C = L/7Z, the operators obey 

{V’i.V’j} = hj, and pj = 

The long wavelength limit of the discrete model cor¬ 
responds to the model in Eq. O while the interaction 
radius 7 Z provides microscopically an ultraviolet cutoff 
in the continuum regime. For 7V-particle states of the 
lattice model we additionally impose the constraint of 
low particle density, N/C < 1, to stay within the con¬ 
ducting regime; a large occupancy N ~ C might lead to 
Wigner crystal physics at sufficiently strong interactions 
that would localise the system. This procedure is ana¬ 
logous to the point splitting regularisation technique^ 
which is usually introduced within the framework of the 
Luttinger liquid mode in the linear regime. 


III. SPECTRAL PROPERTIES 

The model in Eq. ([2]) can be diagonalised via the Bethe 
ansatz approach which is based on the observation that 
the eigenstates are superpositions of plain waves. This 
method is also called coordinate Bethe ansatz^ The ei¬ 
genstates, following Ref. El can be parameterised with 
sets of N quasimomenta kj , 

* = y eiE ‘ kp,jl+i E,<i ' VPi,Pi ' • • • yk i vac > • 

V,ji<—<j N 

(3) 

Their corresponding eigenenergies, H T = FT, are E = 
22^=1 (1 — cos kj) /m. Here |vac) is the vacuum state, the 
scattering phases are fixed by the two-body scattering 
problem, 


+ 2 Vn' 


e^+kv) + 1 + 2 mUe ikl 
e %{ki+k lt ) + l + 2 mUe lk i' 


(4) 


Generally, the system of equations in Eq. (|5j) has to be 
solved numerically to obtain the full spectral structure of 
the observables. However, in the long-wavelength regime 
the solutions can be evaluated explicitly. 

In this limit the scattering phases in Eq. (J^J) 
are linear functions of quasimomenta, 2+/// = 

(ki — kv) / ^1 + (m£/) -1 ^ + 7r, which makes the non¬ 
linear system of Bethe ansatz equations in Eq. (|5j) a 
linear system^ Then, solving the linear system for C 1 
via the matrix perturbation theory up to the first sub¬ 
leading order in 1/C we obtain 


kj 


27 Tlj 
_ mUN 
mU-\-1 


rnU 

mU + 1 


E 




2ttIi 

\”2 

_ mUN \ 

mU +1 J 


(6) 


Note that this calculation is valid for any interaction 
strength at low densities. The corresponding eigenenergy 
and total momentum (protected by the translational in¬ 
variance of the system) are 


£ = E 

j 


2m 


(7) 


and P = Yjj 

The spectrum of the many body states is governed by 
the first term in Eq. ©• Reduction of the quantisation 
length in the denominator of the first term in Eq. ([6]) 
is an exclusion volume taken by the sum of interaction 
radii of all particles. Thus all A-particle eigenstates at 
an arbitrary interaction strength are given straightfor¬ 
wardly by the same sets of integer numbers Ij as the free 
fermions’ states, e.q. the ground state corresponds to 
Ij = -N/2...N/2. 


For example, this result can be used to calculate the 
low energy excitations explicitly that define the input 
parameters of the Luttinger-liquid model, the velocity 
of the sound wave v and of the Luttinger parameter 
K. The first pair of the particle-like excitations, when 
an extra electron is added just above the Fermi energy, 
have J/v+i = 7V/2 + 1 and //v+i = 7V/2 + 2. The dif¬ 
ference in their energies and momenta are — E\ = 


(27 xf N/ 


2 ™(£-^y 


and P 2 —P 1 = 27t/ C. Eval¬ 


uating the discrete derivative, which gives the slope 
of the dispersion around the Fermi energy, as v = 
(. E 2 — E 2 ) / (P 2 — Pi) we obtain 


and 22v i s a sum over permutation of quasimomenta. 
The periodic boundary condition quantises the whole set 
of N quasimomenta simultaneously, 

Ckj — 2 (fji = 2ttI j (5) 

where Ij are sets of non-equal integer numbers. 


» =- 5 , 

(1 NmU \ V C(l+mU)J 

V C{l+mU)) 

(8) 

where vp = 1 tN/ ( mL ) is the Fermi velocity and the re¬ 
lation vK = vp between the Luttinger parameters for 
Galilean invariant systems^ was used. 
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IV. MATRIX ELEMENTS 


A. Algebraic Bethe ansatz 


Now we turn to calculation of matrix elements. But 
first we need to select operators that corresponds to spe¬ 
cific observables. Our interest lies in the dynamical re¬ 
sponse functions that correspond to adding/removing a 
single particle to/from a correlated system and to cre¬ 
ating an electron-hole pair excitation out of the ground 
state of a correlated system. For example, the first type 
of dynamics can be realised in experiments using semi¬ 
conductor nano-structure s 1 h 12 where an electrical cur¬ 
rent, generated by electrons tunnelling into/from the 
nano-structure with their momentum and energy under 
controlled, probes the system. 

The response of the many-body system to a 
single-particle excitation at momentum k and 
energy e is described by spectral function^ 
A(k,s) = —Im [ f dxdte l ^ kx ~ £t ^ G (x, 0, £)] sgn (e — p) /tt. 
Here p is the chemical potential and G(x,x',t) = 
—i (T ( e~ lHt ip (x) e lHt ip^ (V))} is Green function at zero 
temperature. In terms of the eigenstates the spectral 
function reads 


The wave function of the TV-particle eigenstates are 
factorised in the algebraic representation which allows 
the general calculation of various scalar products between 
them. Here we will follow the construction in Ref. E3 
for XXZ spins chains changing basis from 1/2-spins to 
spinless fermions. 

The so called R-matrix acts on a tensor product V 1 GV 2 
space and depends on an auxiliary parameter u, where 
Vi and V 2 are element-element subspaces each of which 
consists of two states |0) • and |1).. It is a solution of 
Yang-Baxter equation R 12 (u\ — U 2 ) R 13 {u\) R 23 (^ 2 ) = 
R 23 (^ 2 ) R 13 {u\) R 12 (ui — U 2 ). For the lattice model in 
Eq. ([2]) the R-matrix reads 

Ri 2 = l-(l-b (u)) (cjci + 4 c 2 ) 

- 26 (u) c\c\c\c 2 + C (u) (c\c 2 + 4 c l) (11) 

where 


A(k,e) = Lj2 | W (0) |0) | 2 5 {e-E f + E 0 ) S (k - P f ) 

f 

+ L J2 I (0|V> (0) I/) | 2 6 (e + E f - Eo) S (k + P f ), (9) 

f 

where Eo is the energy of the ground state |0), and Pf 
and Ef are the momenta and the eigenenergies of the 
eigenstates |/); all eigenstates are assumed normalised. 

Creation of an electron-hole pair out of the cor¬ 
related state at zero temperature at momentum k 
and energy e is described by dynamical structure 
factor— S (fc, s) = f dxdte l ^ kx ~ £i ^ (p (x, t) p (0, 0)), where 
p (x, t) = e~ lHt p (0) e lHt is the density operator evolving 
under the Hamiltonian to time t and the average (...) is 
taken over the ground state. In term of the eigenstates 
the dynamical structure factor reads 

5 (k, e) = L | (f\p (0) |0) | 2 6 (e - Ef) 5 (k - P f ). 

f (10) 
Thus, we will be analysing the expectation values of 
the local operators ^ (0) and p (0). To proceed with this 
calculation we will borrow the result from Ref. U3 and 
El for Heisenberg chains. Our strategy is to perform the 
full calculation for the discrete model in Eq. (|2]) obtaining 
the matrix elements of ipj and pj as determinants. Then 
we will take the long wavelength limit to evaluate the 
form factors for the continuum model explicitly which 
will be the main technical result in the theoretical part 
of this paper. Below we will construct the algebraic form 
of Bethe ansatz, use the Slavnov’s formula^ to express 
the scalar product and the normalisation factors in this 
representation, and finally calculate the matrix elements 
of the local operators. 


b{u) 


sinh (u) 
sinh (u + 2rj) ’ 


sinh (2p) 
sinh (u + 2 rj) 


( 12 ) 


Here p is the interaction parameter, and the tensor 
product space is defined using fermionic basis |0) • and 

\l)j with corresponding fermionic operators j q, cj j = 

Sij that act in these bases as cj |0) • = |1) •. The latter 
will account for anticommuting nature of the lattice fer¬ 
mions on different sites in contrast to the commutation 
relation of the spin operators of a spin chain mode h 22 i 23 
All further calculation are identical to spin chains where 
the anti-commutation relations of the Fermi particles are 
however automatically fulfilled. This approach is more 
convenient than direct mapping of the results for spin 
chains using Jordan-Wigner transformation.— 

A two-states subspace of R-matrix can be identified 
with the two-states fermionic subspace of the lattice site 
j of the model in Eq. ([2]). Then, the quantum version of 
the Lax operator (the so called L-matrix) can be defined 
as Lj = R^j. In the auxiliary subspace £ its matrix and 
operator forms are 


cosh ( u— 77(2 pj —1)) 

cosh(u—77) 

. sinh 2?7ct 

- % ^ 

cosh(u—77) 

= A 7 ^1 — + D^c^c^. (13) 

Here the top left element element of the matrix is a trans¬ 
ition between |0)^ and (0^ states of the auxiliary sub¬ 
space, Cj and pj are the fermionic operators of the lattice 
model in Eq. (j2j), and A- 7 , FF , C\D^ label the matrix ele¬ 
ments of Lj . The prefactor in front of Lj was chosen such 
that for u = Z7r/2 — p it becomes a permutation matrix 
and for p = 0 the L-operator is diagonal. 



. sinh 277c ■ 

— ^_ 1 _ 

cosh(u—77) 

cosh(u+77(2pj —1)) 
cosh(u—77) 
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By construction the L-operator satisfies algebra gener¬ 
ated by Yang-Baxter equation, 

R (u — v) (. Lj (u) (g) Lj (v)) = ( Lj (v) <g> Lj ( u )) R(u — v ). 

(14) 

The entries give commutation relations between the mat¬ 
rix elements of L-matrix. Here we write down three of 
them that will be used later, 




c(u — v) 
b(u — v) 


A3 nj — 

AuCv ~ b(v-u) 


ci A l - C ALZA 

CvAu b(v-u) 


D3 ^—W^) ciDi + W^) 


DiAi), 

(15) 

\ciAi, 

(16) 

-fiDl 

(17) 


A general state of N particles’-Bethe state-is construc¬ 
ted by applying the operator C (u) N times with different 
values of the auxiliary variable Uj , 

N 

* = IHK)I°>, (21) 

i=i 

where a set of N values Uj corresponds to N quasimo¬ 
menta kj in coordinate representation of Bethe states in 
Eq. ([3]). The state in Eq. (f2Tj) with an arbitrary set of Uj 
is not an eigenstate of the transfer matrix r. For instance, 
it can be seen by commuting the operators A and D from 
left to right through all operators C ( Uj ) that generates 
many different states according to the commutation re¬ 
lations in Eqs. (116I17D . However the contribution of all 
of the states that are non-degenerate with T can made 
to be zero by choosing particular sets of Uj that satisfy 
the following set of equations 


These relations can be also be checked explicitly by direct 
use of the definition in Eq. and the Fermi commut¬ 
ation relations. 

The transition matrix T (u) for a chain with C sites-the 
so called monodromy matrix-can be defined similarly to 
the classical problem as 


£ 

T(u) = ^Lj(u). 
J=i 


(18) 


a ( u j ) _ / -i \iV—1 TT- 

d ( u 3 ) diij b (“j - u *) 


( 22 ) 


(see Appendix A for details). 

Under the substitution of the vacuum eigenvalues 
a(uj ) and d(uj ) of A and D operators from Eq. (l20l) 
and b ( ui — Uj) - which define the commutation relations 
Eqs. (1161171) — from Eq. m this so called eigenvalue 
equation above becomes 


If all single-site L-matrices satisfy Eq. d then the T- 
matrix also satisfies the same Yang-Baxter equation, e.g. 
see proof in Ref. 0. Therefore the matrix elements of 

T = A (^1 — in the 2x2 auxiliary 

space £ obey the same commutation relations in Eqs. d 
ED- The transfer matrix for the whole chain, 

r = strT = A(u) — D (u) , (19) 

is the super trace of T-matrix due to the fermionic defin¬ 
ition of the auxiliary space j 22 i 23 The latter gives a fam¬ 
ily of commuting matrices [r (u ), r (?;)] = 0, which con¬ 
tain all conserved quantities of the problem including the 
Hamiltonian. 

The vacuum state |0)—in the Fock space of the model 
in Eq. m -is an eigenstate of the transfer matrix r. The 
corresponding eigenvalue, r (u) |0) = (a (u) — d(uj) |0), 
is the difference of the eigenvalues of the A and D 
operators which can be obtained directly by use of 
the definitions in Eqs. d ttH- Noting that for 
C = 2 Eq. (H8l) gives A(tz)|0) = a\ (u) <22 (u) |0) and 
T)(0)|0) = d\ (u) d 2 (u) |0), where a\(u) = CL 2 (u) = 
cosh (u-\- rf) / cosh (u — rf) and d\ (u) = ^2 (u) = 1, and 
generalising this observation for arbitrary C one obtains 

£ 

, x cosh (^ + 77 ) , 7/ N 

a (u) =-and d(u) = 1. (20) 

cosh (u — rf) 


cosh (uj — rf) C 
cosh ( Uj + rf) 


N 


{-if - 1 n 


sinh (uj — ui — 2rf) 
sinh (uj — ui + 2rf) 


(23) 

Thus all the sets of Uj that satisfy the above equation give 
eigenstates of the transfer matrix in the representation 
of Eq. (|2T1) with the corresponding eigenvalues r (u) T = 
r (u) T where 


r <•■)=° <“> n n i^h-y 


(24) 

This eigenvalue equation in the algebraic framework is 
the direct analog of the Bethe ansatz equation (j5j) in the 
coordinate representation. Direct mapping between the 
two is done by the substitution of 


cosh ( Uj — rf) 
cosh ( Uj + rf) ’ 


mU = — cosh 2ry, 


(25) 


in Eq. © and by taking its exponential. 

!The original lattice Hamiltonian can be obtained from 
the transfer matrix r{u) that contains all of the con¬ 
served quantities of the problem. Logarithmic derivatives 
of t(u) give the global conservation laws by means of the 
so called trace identities, see Ref. m The linear coeffi¬ 
cient in the Taylor series around the point u = ^ — 77 
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is proportional to the Hamiltonian itself. After restoring 
the correct prefactor the expression reads 

H = - S -^d u ln T (u)\ u=if _ v . (26) 

Substitution of the interaction parameter 77 from Eq. ([25]) 
in terms of the particle-particle interaction constant, £/, 
in to the right hand side of the above relation recovers 
the lattice model in Eq. ©• 

B. Scalar product 

The basic quantity, which calculations of expectation 
values will be based on, is the scalar product of two wave 
functions. A general way of evaluating it is the commut¬ 
ation relations in Eqs. (115ill7D and the vacuum expecta¬ 
tion values of the A and D operators. The result of such 
a calculation simplifies greatly if one of the Bethe states 
is an eigenstate of the transfer matrix r(u ), as was first 
shown by Slavnov^ Then the same result was rederived in 
Ref.GjI an d[Ii using the so-called factorising F- matrix^ 
which is a representation of a Drinfeld twist *21 The latter 
will not be used in this Subsection but it will be needed 
later in calculations of the matrix elements of the local 


operators. 

Let |u> = nf= =1 C ( Uj ) |0) be an eigenstate of the trans¬ 
fer matrix so that N parameters Uj satisfy the Bethe 

equation in Eq. (|23|) . And let (v| = (0| YljLi B i. v j) be 
another Bethe state parametrised by a set of N arbitrary 
values Vj. The scalar product of these two states (v|u) 
can be evaluated by commuting each operator B (vj) 
though the product of C (Uj ) operators using the com¬ 
mutation relation in Eq. which generates the A and 
D operators with all possible values of Uj and Vj. They, 
in turn, have also to be commuted to the right through 
the remaining products of the C (Uj ) operators. Finally 
products of the A and D operators, which act upon the 
vacuum state, just give products of their vacuum eigen¬ 
values a(uj ), d(uj ) and a(vj),d(vj ) according to Eq. (1201) . 
The resulting sums of products can be written, using the 
relation between Uj in Eq. (1231) , in a compact form as a 
determinant of an N x N matrix^ 

Ilf; 1 sinh (vj — uA 

<v|u> = ff- • U net -w -7 det 5 

Yljd smh m - Vi) Uja sinh (uj - m) 

(27) 

where the matrix elements are S a b = d Ua T ( v b ). Under 
substitution of the eigenvalues of the transfer matrix from 
Eq. ([24l) these matrix elements read in explicit form as 


cosh L (v b + rf) sinh (2rf) sinh (Uj ■— v b + 2r]) 

cosh L (v b - rf) sinh 2 (u a ~ v b ) .^ a sinh (uj - v b ) 


(- 1 ) 


N sinh (2rj) 


N 


sinh (v b 


u) n 

Ua >j=l^a 


sinh ( v b — Uj + 2rf) 
sinh (v b — Uj) 


(28) 


For N = 1 the result in Eqs. (l27l l28l) follows directly 
from Eq. m • For arbitrary N the proof is more com¬ 
plicated: it employs the residue formula^ (the function 
(v|u) has first order poles when Vi — >> Uj) and the recur¬ 
rent relation for the scalar product of TV + 1 particles in 
terms of the scalar product of N particles, see also details 
in Ref. HU- 

The normalisation factor of Bethe states can be ob¬ 
tained from Eq. (l27l) by taking the limit v —>> u. The 
first order singularities, {v b — u b ) _1 , in the off-diagonal 
matrix elements Eq. ([28]) are cancelled by zeros in the 
numerator in Eq. ([27l) . The diagonal a = b matrix ele¬ 
ments contain second order singularities (v b — u b )~ 2 for 

I 


v —»• u. However, the numerator also becomes zero when 
v —» u in the leading order. Its expansion up to the first 
subleading order cancels the second order singularity of 
the denominator giving a finite expression for the matrix 
elements in the limit. The normalisation factor is found 
to be 


N 

(u|u) = sinh^ (2r]) 


sinh (Uj — Ui + 2rf) 
sinh (Uj — Ui) 


where the matrix elements are 


det Q, 
(29) 


Qab — 


-c 


sinh 2 r] 


£ 


cosh(w a +?7) cosh(n a —77) 
_sinh 4?7_ 

_ sinh(u{,—u a + 277 ) sinh(u{,— u a — 2 rj) 


sinh 4?7 


j7^a sinh(n a —Uj — 2 r]) sinh(u a — Uj +277) 


i a = b, 
5 a 7 ^ b. 


(30) 


The last formula was originally derived by Gaudin us¬ 
ing quantum mechanical identities in the coordinate rep¬ 
resentation of Bethe ansatzj^i Mapping of the resulting 


expression in Ref. [H to the algebraic representation by 
means of Eq. ([25]) gives directly the result in Eqs. (]27l 
1281) with a different prefactors due to different normal- 
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isation factors in the definitions of the states in Eq. © 
and of the states in Eq. (ED). We will use the algebraic 
form in Eq. ED for the calculation of the local matrix 
elements below. 


C. Expectation values of local operators 

Operators of the algebraic Bet he ansatz in Eqs. m 
□§ are non-local in the basis of the original fermionic 
operators of the lattice model in Eq. (|2j). Thus, the first 
non-trivial problem in calculating the matrix elements of 
the local operators f and pi in the algebraic representa¬ 
tion of Bethe states in Eq. ED is expressing the operators 
of our interest in terms of the non-local A, B,C and D 
operators from Eqs. mm- Alternatively these Bethe 
operators can be expressed in terms of the local operators 
of the lattice model. The latter approach is much more 
complicated since the product of matrices in Eq. m is 
a large sum (exponential in the number of sites in the 
chain) restricting severely the ability to do explicit cal¬ 
culations using the fermionic representation in practice. 

An alternative way was found by constructing the F- 
matrix representation of a Drinfeld twist ^ In the F- 
basis the monodromy matrix in Eq. (fl8l) becomes quasi¬ 
local, i.e. its diagonal elements A and D become dir¬ 
ect products of diagonal matrices on each site over all 
sites of the chain and the off-diagonal B and C are 
single sums over such direct products. Direct calcu¬ 
lations become much easier in this basis. Specifically, 
analysis of A, B, ( 7 , D operators leads to a simple result 
for representing the fjj operator in terms of algebraic 
Bethe ansatz operators, which then is shown to be basis 


independent j 17 d 8 

d = rj_1 (y - c (y - r ') rC ~ j (y - *>) • 

( 31 ) 

Here r (u) = A ( u) — D (u) is the super trace of the mono¬ 
dromy matrix and C (u) is its matrix element. 

The transfer matrices in the right hand side of the 
above equation give only a phase prefactor in the expect¬ 
ation values with respect to the Bethe states in Eq. (ED- 
Let |u) be an eigenstate of the transfer matrix with N 
particles, let |v) be an eigenstate with N + 1 particles, 
and let us consider the case of j = 1. Acting with 
the r c ~ x (in/ 2 — rf) operator on the eigenstates |u) gives 
the eigenvalue Y\j=i cosh^ -1 ( Uj — rj) / cosh^ -1 ( Uj + rf) 
according to Eq. ([24l) . Then, using the mapping to the 
coordinate representation in Eq. ED and the Bethe equa¬ 
tion in the form of Eq. ©, this eigenvalue can be ex¬ 
pressed as exp [iP u (C — 1)] where P u is the total mo¬ 
mentum of the state Uj , a quantum number. Similar 
phase factors for j / 1 are evaluated in an analogous 
way and each of them cancels out under modulus square 
in the form factor in Eq. © making the local form factors 
independent of j in full accord with the translational in¬ 
variance of the system and the observable in Eq. (|9j) . Thus 
we will only calculate the value of |u^. 

Since C — rf) Y[J=i C ( u j ) |0) is also a Bethe state 
1^ — 77,^), though it is not an eigenstate, the expect¬ 
ation value can be calculated using the result for the 

scalar product ^v|f|u^ = (v|^ — 77 , i/j). Substituting 
^ — rj,Uj in Eqs. (J27I28D explicitly one obtains 


( 



(_!,»+■ 


. nf=7 cosh ( v j - V) sinh Ar+1 (2r?) det M 

nJLi c ° sh k-+ v) nf«= 2 sinh ( m j - u i ) sinh 


(32) 


where the matrix elements are 


M a b = 


sinh (u b - v a ) 



sinh (u b — Uj + 2rf) 
sinh (ub — Uj — 2rf) 


N+l 

[ sinh (ub 

j = l/a 


- Vi 


AT+1 


277) + n sinh (ub — Vj + 2rf) 

J=l/a 


(33) 


for b < N + 1, and 


cosh (v a - 77 ) cosh (i; a + 77 ) ’ V ' 

for b = N + 1. Here the Bethe equation in Eq. (|23|) 
was used to express a (vj) /d (vj) in the matrix elements, 
and some factors in the matrix elements and the over¬ 
all prefactor cancel out. This result can be checked by 
numerical evaluation of the sums over spacial variables 
using the coordinate representation in Eq. (|3]) for a small 
number of particles N = 1,2,3 which we have done. 


The determinant results in Eqs. ©E ED G© can be 
checked by numerical evaluation of the sums over spacial 
variables using the coordinate representation in Eq. ©• 
However the latter summation over many coordinates has 
factorial complexity which already limits numerical cal¬ 
culations to a few particles on chains of a few dozens 
sites. The results of the algebraic Bethe ansatz calcu¬ 
lations have a power-law complexity that allows general 
studies, at least numerically, of systems with hundreds 
of particles on arbitrary long chains without making any 
approximations, e.g. the studies of correlation functions 
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in one-dimensional systems in Refs. [25l~[30l 


D. The long wavelength limit 


We now turn to the evaluation of the long wavelength 
limit for matrix elements in the determinant form with 
the aim of calculating the determinants explicitly. The 
resulting expressions will then be used to study physical 
observables. 

Such an analysis is more convenient in the coordinate 
representation. For small kj the non-linear mapping in 
Eq. (|25l) becomes linear, similarly to Bethe equation in 
Eq. ([5]) in this limit. Then, a simple inversion of the 
linear function gives 


Uj = ^\l ^\ kj and 71 = ~ ^ acosh ( mC/ ) • ( 35 ) 


Note that \uj\ and kj are simultaneously much smaller 
than one, while the interaction strength U can be of an 
arbitrary magnitude. 

We start from the expansion of the normalisation 
factor in Eq. (|29|) up to the leading non-vanishing or¬ 
der in the quasimomenta. We first substitute Eq. ([35]) in 
the matrix elements in Eqs. then expand them up 

the leading non-vanishing order in kj <C 1 , and obtain 
the diagonal matrix elements as follows, 


Qaa = 2 C 


mU — 1 
mil + 1 


2 (N — 1 ) mU 
Vm 2 C/ 2 - l" ’ 


(36) 


and 


Qab — 


2 mU 


\Jm 2 U 2 — 1 


(37) 


for a 7 ^ b. The off-diagonal matrix elements are small 
compared to the diagonal entries as Q a b/Qaa ~ l/£ so 
the leading contribution to the determinant is accumu¬ 
lated on the diagonal. Also expanding the prefactor of 
Eq. (|29l) in small kj we obtain the following expression 
for the normalisation in the long wavelength limit, 


(k|k> 


2 n ~ (- 


(-if (1 -mur (C-*™) 

i N(N ~ 1 ) n ijij (kj - h) 


N 


(38) 


where kj are quasimomenta in the coordinate represent¬ 
ation of Bethe ansatz. 

Our primary interest lies in the spectral function which 
contains the local matrix element of operators so here 
we will focus on the determinant result in Eq. ([32l) . Sim¬ 
ilarly to the calculation of the normalisation factor we 
substitute Eq. (|35l) into Eq. (l33l) , which however becomes 
zero in the zeroth order in kj. Expanding it up to linear 


order in kj we obtain 

AT-1 k u - V Ar+1 k v 

M ab = 2 mU (mW - l)^~ ^ ^ J 

(39) 

for b < N + 1, where A P = JA k'j — J2j ^ is the dif¬ 
ference of two conserved quantities, the momenta of two 
states k. u and k v . The matrix elements in Eq. (l34l) are 
already non-zero in the zeroth order in kj giving 

Mab = mU + 1 ’ 

for b = N + 1. Also expanding the prefactor in Eq. 
(|32l) and rearranging the expressions by taking a common 
factor out of the matrix elements we obtain 

<kV (0) |k“) = (~1) N+1 i N2 2 N2+N +i 

(mU -l) N2+ i m N U N V 

lx nN {kf-kY) nff 1 (*?-*?) 

where the entries of the matrix under the determinant, 
V = det Ad, for b < N + 1 are 
A P 4- k v 

Mah = Q and M a ,N+i = 1. (42) 

All matrix elements are of the the same order so the 
determinant in Eq. (|4T|) is a sum of a large number of 
terms unlike the normalisation factor in Eq. (]38|) . Doing 
the summation we find an explicit expression in the form 
of a fraction of two polynomials in quasimomenta of the 
initial and the final states, 

„ , ,, w n,(AP+t?) 

n,,, - Ht) 

N AT+1 

II( fc “- fc b II ( k J- k i)- ( 4 3) 

j<i j<i 

For N = 1 the result above is evaluated straight¬ 
forwardly as a determinant of a 2 x 2 matrix with 
the matrix elements in Eq. (l42l) . For arbitrary N 
we prove it by induction. Using the Laplace develop¬ 
ment on the N + 1 row, the determinant for N + 1 
particles can be expressed as a sum of minors given, 
in turn, by determinants for N particles, Xhv+i = 

Yla=i (-l) JV+1+ “ A^ a ,JV+iminor a) jv+i, which - let us as- 
sume for purposes of the inductive method - are given 
by Eq. (1451) 


minor a5 Ar+i = (-1) 


iv+i n; =1 (AP+fc“) 

t AT,AT+2 


I riV ,iV -\- 2 , 

IU=1J=1 #a \ K j ) 

N AT+2 

n (kJ-K) n (kj~k«) (44) 

j<i j<i^a 


Here Ad a? jv+i are given by the matrix elements in Eq. 
62), N quasimomenta k'j are labeled by j = 1... A, and 
N -f 1 quasimomenta k v - are labeled by j = 1... TV + 2 
with a th elements excluded. Under taking common factor 
in front of the sum, the determinant for N + 1 particles 
becomes 
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V 


AT+l 


= (“I) 


AT+2 


n.r (ap 


^) n 


AT+l 

j<i 


(*?-*?) riE 2 (q-kv) 


11,, (kv - /.•") 


1 


" +2 (ap+ k ) njL + i„ (*? - ^+0 nli (*? - 


AP + A: 


E 


A^+l a=l 


nf =1 (*y - ** +1 ) n;= + E (*? - *s) 


(45) 


The sum in the above expression gives, by direct calcu¬ 
lation, YlH=i • • * = AP + k u N+1 which makes the whole 
second line unity. The determinant is equal to the first 
line of Eq. (H5l) which is also equal to the result in Eq. 
fi3l) for N + 1 particles. Thus we obtained the same 
result for TV + 1 particles starting from Eq. m f° r a 
particles. Hence, it is proved by induction. 

Finally, the form factor in the Eq. © is the modulus 
squared of Eq. dnji. Normalising the initial and the final 
state wave functions using Eq. CU) as (0) |o)| = 

(0)|k°)| 2 (k/|kq 1 (k°|k°) 1 we obtain 


(o)|o )| 2 


z nffe 0 -^/) 2 
' c n,"i w+1 

N N+l 2 

n«-‘.°) 2 n(‘/-‘0 • 

i<j i<j 


(46) 


where Z = mU / (mU + 1) / (£ — NmU / (1 + mU)), fcj 
and kj are the quasimomenta of the eigenstate |/) and 
the ground state |0), and Pq = 0 for the ground state. 

Calculation of (/|p(0)|0) is done in a similar way 
by expressing the local density operator pi, within the 
framework of the lattice model, in terms of the algebraic 
Bethe ansatz operators A, P, C, D and, then, by using 
the Slavnov formula. Details are given in appendix B. In 
the long wavelength limit we obtain 


y2N- 

I(/Ip(o) |o)| 2 = 


td2N 


C 2 t~tN,N 
I T,7 


(*/-*?) 

N 2 

n^-^) 2 n(d^0 > ( 4? ) 


i<j 


i<j 


where the final states |/) have the same number of excit¬ 
ations N as the ground state |0), unlike in Eq. (pf6lh and 
Po = 0 for the ground state as in Eq. flU). 

These form factors in Eqs. mm together with the 
solution of Bethe equations in Eq. © is the main tech¬ 
nical result in the theory part of our work. We will 
analyse its physical consequences in the next two Sub¬ 
sections. The similarity between these two expressions 
means that the hierarchy of modes we will identify below 
is a general feature of one and two body operators. 


—(J>00<K>(0<>0-0-C|>— gs) 0 

—kp kp 



Figure 2. Configurations of quasimomenta that solve the 
Bethe equations in Eqs. (mu for the spinless fermions model 
in Eq. ([2j) : gs) the ground state, a) excitations that form the 
a-level of the hierarchy, and b) excitations that form the b- 
level of the hierarchy. 


E. Hierarchy of modes 


The results in Eqs. mm have one or more singular¬ 
ities when one or more quasimomenta of an excited state 
coincide with a quasimomentum of the ground state, 
kj = kj. Both results have a multiplicand Z 2N ~ C~ 2N 
that becomes virtually zero in the thermodynamic limit, 
in which C oo. Thus the product of these two opposite 
factors produces an uncertainty in the limiting behaviour 
(of the 0 x oo type) that has to be resolved. Since we are 
specifically interested in a transport experiment in this 
paper, in which the spectral function is measured, we 
will mainly focus on solving the uncertainty problem for 
result in Eq. (®. 

The maximum number of singularities is N in the ex¬ 


treme case, when quasimomenta kj of an excited state 
coincide with all of the N quasimomenta of the ground 
state kj given in Fig.[2^s. The excited states of this kind 
are given in Fig. [2^. The divergences in the denominator 
of Eq. (H6D occur only in the leading order—the first term 
in Eq. —but the subleading order—the second term in 
Eq. ([6])— already provides a self-consistent cutoff within 
the theory. The interaction shift of the quasimomenta 
at subleading order does not cancel for the extra added 
particle in the excited state, making the factors in the 
denominator of Eq. (H6l) 


k i 


mU kf, +I - k° 


k u = 

J mUPiC 


mUN ’ 
mU +1 


(48) 
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where in the r.h.s. only the first term from Eq. (| 6 ]) is 
relevant for k^ N+1 and k®. The numerator for the states 
in Fig. [2^ becomes 

k j ~ P f= k j ~ k N+v ( 4 9) 

Substitution of Eqs. (|48l l49|) in Eq. (|46l) for one particle, 
say for j = N, cancels one factor Z 2 ~ C~ 2 and the 
other part of the product for i ^ j in the denominator of 
first line of Eq. (|46l) cancels partially the products in the 
second line of Eq. fl®. The expression for the remaining 
N — 1 particles is the same as Eq. (|46l) but the numbers 
of terms in the products are reduced by one, N —»• N — 1 , 
giving 



2 Iff" 1 (*? - P/)” 

N—l N 2 

n ( fc j°- k i ) 2 n (d ~ k () ■ ( 5 °) 

i<j i<j 

Repeating the procedure N—l times we cancel the re¬ 

maining Z 2N ~ 2 factor completely (with the rest of other 
terms) and obtain 

£|</|Vh0)|0)| 2 = l. (51) 

Corrections to this result originate from higher sublead¬ 
ing orders in the solutions to Bethe equations in Eq. (j 6 j) 
and are of the order of O (£ _1 ). This becomes much 
smaller than one at leading order of Eq. (|5ll) in the ther¬ 
modynamic limit. 

Substitution of Eq. in Eq. (j9]) gives the value of the 
spectral function A(k,e(k)) = 1. The energies and the 
momenta of the excitations in Fig. [ 2 ji form a single line 
on the spectral plane, like a single particle with disper¬ 
sion 5 (k) = k 2 / ( 2 m*), where the effective mass is renor¬ 
malised by the Luttinger parameter K , m* = mK 41 
Note that, since we still resolve individual levels here, the 
delta functions in the definition of the spectral function in 
Eq. ([9|) become discrete Kronecker deltas. Thus, A(k,e) 
at each discrete point &, 5 describes the probability of 
adding (removing) a particle, which is non-negative and 
is bound by one from above, instead of the probability 
density as in the continuum case. Dimensional analysis 
makes this distinction clear immediately. 

The excitations that have one singularity less (N — 
1 in total) can be visualised systematically as an extra 
electron-hole pair created in addition to adding an extra 
particle, see Figj 2 j 3 . Staring from Eq. (H 6 l) and using the 
same procedure as before Eq. (1501) but N — l instead of 
N times we obtain 



|</IV + (o) |o )| 2 


Z 2 (fc/-fcf) 2 (fcO_ P/ ) 2 

-4-2> ( 52 ) 

(k( - k i ) (d - k i ) 


Figure 3. The spectral function for interacting spinless fer¬ 
mions in the region — /cf < k < kp (kp < k < 3kp) labelled 
by 0(1). The grey areas mark non-zero values. The green 
and the blue lines are modes of the hierarchy labelled as fol¬ 
lows: p (h) shows the particle(hole) sector, /cf is the Fermi 
momentum, a, 6, c respectively identify the level in the hier¬ 
archy in powers 0, 1, 2 of 7 Z 2 /L 2 , and (r, l) specifies the origin 
in the range—modes on the edge have no such label. 


where k{ , k£ and fcj are positions of two particles and 
one hole in Fig. Et>. Substitution of Eq. (15^1) in Eq. (|9]) 
gives values of the spectral function A(k,s) ~ C~ 2 that 
smaller than the values for the excitations in Fig. [ 2 ji (in 
Eq. ©) by a factor of C~ 2 . For two singularities less 
(N — 2 in total) we find A (fc, e) ~ C~ 4 and so on. 

This emerging structure separates the plethora of 
many-body excitations into a hierarchy according to the 
remaining powers of C~ 2 in their respective form factors. 
We label the levels of the hierarchy as a, 6 , c reflecting 
the factors C~ 2n with n = 0, 1, 2 . While the lead¬ 
ing a-excitations form a discrete single-particle-like dis¬ 
persion, see hOa and pla(l) lines in Fig. [3j the spectral 
properties of the subleading excitations are described by 
a more complicated continuum of states on the energy- 
momentum plane. We will explore the 6 -excitations be¬ 
low. 

All simple modes, formed by single particle- and hole¬ 
like excitations of the type in Fig. [2 Jd, in the range 
—kp < k < 3kp are presented in Fig. [3j We use the fol¬ 
lowing naming scheme: p (h) indicates the particle (hole) 
sector, 0 ( 1 ) encodes the range of momenta — kp < k < kp 
(fcp < k < 3fcp), &, 6 , c reflect the terms C~ 2n with 
n = 0, 1, 2 . The suffix (r) or (r) marks a particle-like 
mode, e.g. the states with in Fig.[ 2]3 with k{ = —kp — 7 , 
k^ = fcp, and fcp > Pf = —2kp + k^ > —kp forms the 
mode p0b(l). Hole-like modes have no suffixes, e.g. the 
states in Fig. (2 Jd with k( = —kp — 7 , k£ = kp + 7 , and 
—kp < Pf = — k^ < kp form the mode pOb. Simple modes 
formed by excitations of lower levels of the hierarchy are 
obtained by translation of the 6 -modes constructed in 
this paragraph by integer numbers of ± 2 kp- A couple 
of simple modes formed by c-excitations are presented 
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x = 0 

X — 1 

pxa 

- 

1 

hxa 

1 

- 

pxb 

16Z 2 k^k 2 

TzvF-F+ItJ 2 

(ft; 2 —O f +7) 2 )^ 

(A; —fc F +7) 2 0 —fc F +27 ) 2 

pxb ( l ) 

4Z 2 O f +/c) 2 

- 

pxb (r) 

AZ 2 (k F -k) 2 

k 2 

- 

hxb 


4Z73fc F —fc —7)7fc F +fc7 

fc£(fc-fc F +7) 2 

hxb ( l ) 

4ZV 

7^ 

Zj rb -p /v 

0 + fc F +27 ) 2 

((fc+ 7 ) 2 -fc F y 

hxb (r) 

4Z 2 7 2 

(k — kp—2'y) 2 

- 


This immediately moves such states one step down the 
hierarchy to 6 -excitations. Let us consider the spectral 
function as a function of energy only making a cut along a 
line of constant k. The energies of the electron-hole pairs 
themselves are regularly spaced around the Fermi energy 
with slope up. However, the degeneracy of the many- 
body excitations due to the spectral linearity makes the 
level spacings non-equidistant. We smooth this irregu¬ 
larity using an averaging of the spectral function over 
energy, 

— f ^ de 

A(k,s) = / — A(k,e-\-e) (55) 

J- e A- e o 


Table I. Spectral weights along the a- and the 6 -modes for 
—kp < k < kp (&f < k < 3/cf) labeled by x = 0 ( 1 ). 
Terminology is the same as in Fig. [3] 7 = 27r/£ and Z — 
mU/ (mU + !)/(£- NmU/ (1 + mU)). 


on Fig. [3] They have the same naming scheme as the 
6-modes. 

Now we evaluate values of the spectral function along 
all simple 6-modes in the range — fcp < k < 3fcp- Let 
us start from the pOb mode, see Fig. [3j Along this 
mode the spectral function is a bijective function of fc, 
A(k,£ p ob (k)) where e p ob (k) = fcp/ ( mK ) — k 2 / (2 mK). 
The states that form it belong to 6-excitations in Fig. [ 2)3 
with k{ = — fcp — 7 , fcjf = + 7 , and k = Pf = —k®. 

Substituting this parameterisation in Eq. (f52l) we obtain 

. /7 /7N . 1 6Z 2 k??k 2 

A ( ft, £ p0b (ft)) = — - -2 • (53) 

(fc2-(fc F+7 ) 2 j 

The spectral function along all other 6-modes in Fig. [3] is 
calculated in the very same way and the results (together 
with a-modes) are summarised in Table HI 

The amplitude of the subleading 6-excitations does not 
vanish in the thermodynamic limit, though it is propor¬ 
tional to 1 /C. The limit involves both C —>> 00 and the 
particle number TV —>> 00 but keeps the density N/C fi¬ 
nite. The spectral weights of the subleading modes p06, 
ft 16, and hlb (r) from Table [T] are proportional to the 
density squared for some values for fc, e.g. the modes pDb 
at k = fcp gives 


A (fcp, £pOb (&f)) 


/ mU \ 2 TV 2 

V 1 + m U) (r NmU A 

1 +mU) 


see Table U for other modes, and are apparent in the 
infinite system. 

Assessing further the continuum of 6 -excitations we 
consider the spectral function and how it evolves as one 
moves away slightly from the of the strongest a-mode. 
Just a single step of a single quantum of energy away from 
the a-mode requires the addition of an electron-hole pair 
on top of the configuration of quasimomenta in Fig. [2ji. 


where eo is a small energy scale. 

Then, using the parametrisation of 6 -excitations in Fig. 
E3> in the vicinity of the principal parabola, we linearise 
the energies of the extra electron-hole pairs around the 
Fermi energy and of the particle around its original po¬ 
sition. We then substitute the resulting expressions for 
k( 2 and k® in terms of the energy E from Eq. © in 
Eq. (|52l) . similar to our procedure of obtaining Eq. ([53]). 
Finally, we use the averaging rule in Eq. (|55 | ) and obtain 


A(k,e) 


A(k,e) 


Zpkp (3fc 2 + fc 2 ) 9 (gfrQa (fc) - e) 


m^K (< s h0a (k) - e) 

for — < k < fcp, 

Z 2 (fc + sgn (e - £pi a (;) (fc)) fcF ) 3 
m^K\e- e pla{l) (ft) | 


(56) 


(57) 


for fcp < k < 3kp 


where 7 = 2tt/C and £hoa (ft) = £ p i o (0 (ft) = ft 2 / (2 mK) 
is the parabolic dispersion of the a-mode. 

The result in Eqs. ( 1551157)1 can be interpreted as the 
line-shape of the a-mode. However, it has an unusual 
form—namely that of a divergent power-law. The diver¬ 
gence at the parabola is cut-off by the lattice spacing 
recovering A(k,£hOa,pia(k)) = 1 from Eq. ([5l]) . In Eq. 
(I57|) the line-shape is asymmetric due to different pre¬ 
factors ( k dz /q?) 3 above and below the line. In Eq. (|56l) 
the higher energy part (e > £ho a (k)) is absent due to the 
absence of the excitation in this region, forbidden by the 
kinematic constraint. 

Not every simple mode marks a distinct feature. The 
states at least on one side, above or below the mode in 
energy, have to belong to a different level of the hierarchy 
than the mode itself, which results in a divergence or in 
a jump of the spectral function in the continuum of ex¬ 
citations. Otherwise the spectral function is continuous 
across all of the modes that belong to the same level of 
the hierarchy as the excitations around them. The a- 
modes are distinct since excitations around them belong 
to a different 6 -level. All modes on the spectral edges, 
p06, plb, hlb , and so on, are distinct since on one side 
there are no excitations (due to the kinematic constraint) 
and on the other side there is a finite density of states 
resulting in a jump of the spectral function. 
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An example of an observable subleading mode in the 
continuum is hOb (r). On the higher energy side of this 
mode, the excitations are described by the same type of 
states in Fig. [2 Jd but on the lower energy side creation 
of an additional electron-hole pair in the quasimomenta 
results in states that have two non-cancelled singularities 
in Eq. (1351) . which lowers their corresponding level of the 
hierarchy to c from b. This, in turn, results in an observ¬ 
able feature in the spectral function at the position of the 
hOb (r) mode. On the other hand, the pOb (r) and hlb (/) 
modes in continuum are not detectable since excitations 
on both sides around them belong to the same 6-level of 
the hierarchy. Observability of all other modes can be 
easily assessed in the same way by considering their cor¬ 
responding states in the form of Fig. [2^ and Fig. [2 ]d and 
excitations around the modes. 

The structure of the matrix element in Eq. <137)1 is quite 
similar to the matrix element in Eq. fl6l) . Thus the dy¬ 
namical structure factor exhibits the same hierarchy of 
excitations (and modes formed by them) as the spec¬ 
tral function analysed in detail in this Subsection. The 
strongest excitations correspond to only a single electron 
hole-pair, the first subleading level corresponds to two 
electron-hole pairs, and so on. 

Also, a similar hierarchy of excitations was observed in 
numerical studies of spin chains, e.g. Refs. [25l [27l . ItO [33l . 
There, it was found that only a small number of electron- 
hole pairs are sufficient to saturate the sum rules for the 
dynamics response functions. For example, integration 
of Eq. (|§]) over the energy and momentum, 


A, 



(Ej +1 - E j )2mK/k F /'y 


Figure 4. Distributions of level spacings in the vicinity 
(max (Ef — € p ob (k)) /ef — 1/100) of the particle mode pOb 
accumulated along energy axis for the values of momenta a) 
k = 0 and b) k = 0.4355/c F ; N = 2 • 10 3 and L = 2 • 10 5 . 


behaviour A(k,s) ~ \e — £ e dge(&)| see Refs. 0IM- 
h7 MM For spinless fermions the exponent of the spec- 
tral function is given by^L 



(59) 


dedkA (fc, e) = C — N, (58) 

gives the number of empty sites. If a sum over only a 
small number of electron-hole pairs in the intermediate 
state / in Eq. ([9 ]) is sufficient to fulfil this rule in Eq. ([58] ) 
then a few electron-hole pairs already account for major 
part of all of the spectral density and the states with 
more electron-hole pair have vanishing spectral weights, 
as in the hierarchy of modes established in this work. 
Our analytic work demonstrates how this can arise in a 
Bethe Ansatz solution, though the numerical studies of 
spin chains were done at large fillings (C ~ N ), for which 
our result in Eq. (|46l) is not directly applicable. 


F. Spectral edge modes 

In this Subsection we consider another important role 
played by the continuum of eigenstates, namely how they 
form the spectral function close to the spectral edges. 
These edges separate regions where there are excitations 
from regions where there are none, see borders between 
white and grey regions in Fig. [3j The recently proposed 
model of a mobile impurity^— gives a field-theoretical 
description of the dynamic response functions around the 
spectral edges predicting a general (divergent) power-law 


for both the particle (pOb) and hole edges (hOa), where 
Eq. (|5J) gives the Luttinger parameter K in terms of 
the microscopic parameter of the model in Eq. ©• Here 
we will compare the field-theoretical result in Eq. (|59l) 
with the microscopic calculation in Eqs. ® S6D. We 
find agreement in many cases, but interestingly we also 
find some cases where the mobile impurity results are 
not consistent with the analytic solution, suggesting this 
field-theoretical approach is not the complete story. 

The hole edge is an a-mode, hOa, whereas the con¬ 
tinuum around it is dominated by 6-excitations. The 
spectral function formed by these 6-excitations has 
already been calculated in Eq. HMD giving the power-law 
behaviour with the exponent a = 1. Note that for spin¬ 
less fermions the Luttinger parameters have only small 
deviations from K = 1 for arbitrary magnitude of the 
short-range interactions. This makes the result in Eq. 
(I59| ) a = 1 for all values of U ; small deviations (which 
are 67-dependent) make a < 1 and require subleading 
terms in 1/C in the Bethe ansatz calculation in Eqs. (J6J 
USD. Thus the result of the microscopic calculation coin¬ 
cides with the prediction of the mobile-impurity model 
in Eq. ([59] ) for the hole edge. 

The particle edge is a 6-mode, p06, and the excitations 
around it belong to the same 6-level of the hierarchy as 
the mode itself. Parameterising the 6-excitations in this 
region of the continuum as in Fig. [2}d and using the the 
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averaging procedure in Eq. (l55l) we obtain, repeating the 
same steps as before, Eqs. mm, 

A (e) ~ (e - s p ob ) 3 (60) 

for k « 0 to 

A (e) ~ const — (e — e p ob) (61) 

for k « fcF, where e p o b(k) = k\f ( mnK ) — k 2 / (2 mK). 
This is a new power-law behaviour characterised by an 
exponent a changing essentially with k from a = — 3 
for = 0 to a = — 1 for k ~ zb/cp and is different 
from the predictions of the mobile-impurity model in 
Eq. (1591) . Here we observe that the phenomenological 
model in Refs. [334361 is correct only for the a-mode spec¬ 
tral edge but higher-order edges would require a different 
field-theoretical description. 

On a more detailed level, difference between the 
particle and the hole edges manifests itself in dif¬ 
ferent statistics of level spacings around the edges. 
Evaluation of the density of states, v(k,e) = 
8 ~ Ef) 8 (& — -P/), is performed using Ef from Eq. 

© for a fixed momentum k. For 6-excitations in Fig. [2)3 
we obtain the same results, 

v{k,e) ~ |e - £ p0b (hOa) (&)|. (62) 

in the vicinity of both the particle pOb and the hole hOa 
edges. However the statistics of the level spacings 

V(s,k) = J2Hs- (Ef +1 - Ef)) S (k - Pf) , (63) 

f 

where Ef are assumed sorted by their values, is different 
in the two regions. For the hole edge the energy levels 
are spaced regularly and are governed by the slope of dis¬ 
persion at the Fermi energy ~ v . This gives a bimodal 
V (s, k) with a sharp peak at s = 0 (due to many-body de¬ 
generacy of almost linear spectrum at Fp) and at s ~ vj. 
For the particle edge the statistics of the level spacings 
varies from having a regular level spacing (for k commen¬ 
surate with kp in Fig. [4^) to an irregular distribution (for 
incommensurate k in Fig. HJd). 

The change in the characteristics of the underlying 
statistics is another microscopic difference between the 
particle (6-mode) and the hole (a-mode) edges that sig¬ 
nals a difference in underlying physics for the particles 
and for the holes spectral edges beyond the low energy 
region. 


V. LOCAL DENSITY OF STATES 

Now we turn to another macroscopic observable, the 
local density of states (LDOS), which describes the prob¬ 
ability of tunnelling a particle in or out of the wire at a 
given position in space and at a given energy. Since the 
model in Eq. © is translationally invariant the LDOS 



Figure 5. The local density of states for interacting spin¬ 
less fermions: the red and the green lines show the contri¬ 
bution of a- and 6-excitations and the blue line indicates the 
Luttinger-liquid regime. Inset is a log-log plot around the 
Fermi energy: the points are numerical data for N = 71, 
L — 700, mV = 6 giving K = 0.843, and the dashed line is 
n (e) = const \e — p\ ( K+K )/ 2_1 . 


depends only on a single variable-energy, making it a 
more convenient quantity to study qualitatively how the 
physical properties change from low to high energies. In 
this section we will show how the power-law result of the 
Tomonaga-Luttinger model^ at low energy crosses over 
into the hierarchy of modes dominated behaviour at high 
energy. 

The probability of local tunnelling at en¬ 
ergy e and at position x is described byi£ 
n(x,s) = — Im[ f dte~ l£t G(x, x, £)] sgn(£ — p)/tt 
where p is the chemical potential and G(x,x',t) = 
—i (T (e~ lHt f> (x) e lHt ^ (#'))) is the two point cor¬ 
relation function at zero temperature. In terms of 
eigenmodes it reads 


n (e) = C Y [ I (/ V (0) |0) | 2 5 (e-E f + E 0 ) 
f 

+ 1(01^ (0) |/)| 2 S (s + Ef — Eq) 


(64) 


where the coordinate dependence drops out explicitly, the 
eigenenergy Ef have already been calculated in Eq.© 
and the matrix elements |(0|^ (0) |/)| 2 is in Eq. (pf6l) . 
Note that the definition in Eq. (IM1) is connected to the 
definition of the spectral function in Eq. m via 

n{e) = J dkA(k,s). (65) 

The leading contribution for 5 > 0 comes from a- 
excitations. Substituting the matrix element for the a- 
excitations from Eq. ([5ll) we sum over the single-particle- 
like excitations (with e = k 2 / (2 mK)) that form the 
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mode and obtain 


, x / 2 mK . , v , v 

n(e) = y-0 (e). (66) 

This result gives the same 1/v^ functional dependence 
-see red line in Fig.[5]-and the same 1 /^ van Hove singu¬ 
larity at the bottom of the band £ = 0 as the free-particle 
model. 

For 5 < 0 the leading contribution to n (e) comes from 
6 -excitations. Instead of performing a summation in Eq. 
(ED over every 6-excitation in this region we use an inter¬ 
mediate result in Eq. (|56l) where the matrix elements of 
6 -excitations in Eq. (IS2J) are already smoothed over many 
eigenstates and the relation in Eq. (f65l) . Evaluating the 
integral over k in Eq. ([65lh after the substitution of Eq. 
(j52l) into it, for e < 0 we obtain 


n(e) 


2 z^4 

7 fiK 


9 (s) 


X 



(67) 


However, away from the linear region the particle-hole 
symmetry of the Tomonaga-Luttinger model is broken 
by the finite curvature of the dispersion and only the 
leading 1/C order in Eq. ([46]) is sufficient to account for 
the main contribution there. 

The general picture emerging in Fig. [5] is a power-law 
crossover between different energy scales. At low ener¬ 
gies (blue region in Fig. ED Eq. (1461) breaks down and 
the collective modes of the Tomonaga-Luttinger model 
are a better description of the system. At high energies 
(the red and the green regions in Fig. ED the hierarchy of 
modes, which directly follows from Eqs. (nail), becomes 
the dominant physical picture. For spinless fermions the 
extent of the crossover region is large due to only small 
deviations from K = 1 for arbitrary short-range interac¬ 
tions. For very small exponents \[K + if -1 ) /2 — l] < 1 
the power-law in Eq. deviates significantly from 1 
only in an extremely narrow region around e — hav¬ 
ing a large window of energies where it overlaps with the 
a-mode result in Eq. (f66l) . 


VI. EXPERIMENTS ON SPINFUL FERMIONS 


There is a finite probability to find a particle below 
the bottom of the conduction band-green line in Fig. 
[5]-which is allowed only due to interactions between many 
particles. The factor Z is proportional to the interaction 
strength V-see Eq. flU) -making n (e) = 0 for e < 0 in 
the free particle limit of V = 0. At the bottom of the 
band below the e = 0 point in Fig. El the result in Eq. 
([67]) contains another Van Hove singularity 




27rZ 2 fc| 

7-^vTFl 


( 68 ) 


which also disappears when V = 0. The appearance of 
the identical exponent as in Eq. ([66]) exponent 1 / \/\s\ 
seems to coincidental. 

Around the Fermi energy (the point 5 = fi in Fig. 
© the Tomonaga-Luttinger model predicts a power-law 
suppression of LDOS, 


n (e) rsj \s — /i\( K+K 1 , (69) 

e.g. see the book in Ref. @. However the result f or the a- 
mode in Eq. ([66]) is finite at this point, n (/i) = ^2mK/ /x. 
In order to resolve the apparent discrepancy we evaluate 
n (e) numerically around the Fermi energy using the de¬ 
terminant representation of the form factors Eqs. (|32J[33D 
instead of Eq. (|46l) , which accounts for all orders in 1/C, 
and indeed find a suppression of LDOS around e = /i, 
see blue line in the inset in Fig. [5j This signals that 
the leading-order expansion in the C |(/|V^ (0) |0)| = 1 
result is insufficient at low energies. Very close to the 
Fermi point (the linear region of the single-particle dis¬ 
persion) all 1/C orders of the Bethe ansatz calculation 
are needed to reproduce the result of the Tomonaga- 
Luttinger model, see dashed lines in the inset in Fig. [5] 


So far, in this paper, we have established the theoret¬ 
ical framework for expecting a hierarchy of modes in a in¬ 
teracting system at high energy. Now we turn to a meas¬ 
urement of tunnelling of electrons in a one-dimensional 
(ID) nanostructure, which gives experimental evidence 
for the existence of the hierarchy. Electrons have spin 
1 /2, which does not correspond directly to the model of 
spinless fermions in Eq. ©, and there is currently no 
known method for calculating the necessary form factors 
for spinful fermions. However, the general picture that 
emerges from the experiment is qualitatively the same as 
our result in the theory part of this paper. 

The design of our device^ is based on a high-mobility 
GaAs-AlGaAs double-quantum-well structure (blue and 
yellow layers in Fig. ED, with electron densities around 3 
and 2 x 10 15 m -2 in the upper and lower layers, respect¬ 
ively, before application of gate voltages. Electrons in 
the upper layer are confined to a ID geometry (‘wires’) 
in the x-direction by applying a negative voltage to split 
‘finger’ gates on the surface (gold layer in Fig. [6j. 

Electrons underneath the gates are completely de¬ 
pleted, but electrons below the gap between gates are 
squeezed into a narrow ID wire. The extremely regu¬ 
lar wires are arranged in an array containing ~ 600 of 
them to boost the signal. The small lithographic width 
of the wires, ^ 0.18 /im, provides a large energy spacing 
between the first and second ID subbands defined by 
spatial modes perpendicular to the wires (^3 — 5 meV, 
probably somewhat smaller than for overgrown wiresii) . 
This allows a wide energy window for electronic excit¬ 
ations in the single ID subband that covers a range of 
a few chemical potentials of the ID system. The lower 
2 DEG (blue in Fig. ED, is separated from the wires by a 
d = 14 nm tunnel barrier. The wafer is doped with Si 
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Figure 6. Schematic of the device made out of a double¬ 
well heterostructure. The dark blue and cream layers are the 
lower and upper quantum wells, respectively. The lower layer 
hosts the two-dimensional electron gas (2DEG). The wires 
are defined in the upper layer by gating. The gold top layer 
represents metallic gates deposited on the surface of the semi¬ 
conductor heterostructure. The array of parallel ‘finger’ gates 
defines the ID wires in the upper well. The white lines rep¬ 
resent air bridges joining the finger gates together. Current 
is injected from the ohmic contact on the right solely into the 
upper well through the constriction at top right in the dia¬ 
gram. The constriction is formed and pinched off by a split 
pair of gates, and charge is induced again in the upper well in 
the constriction by a gate in the centre of the channel. The 
current then flows into the ID wires via the narrow, nomin¬ 
ally 2D, regions shown in light blue. Tunneling to the 2DEG 
below is possible, and this gives a small ‘parasitic’ current in 
parallel with the tunnel current from the ID wires. To tune 
this tunneling off resonance in the regions of interest by chan¬ 
ging the density, a ‘p’ gate is placed above the ‘p’ regions and 
a voltage Vp is applied. Current is prevented from flowing 
from the upper well into the left-hand ohmic contact by a 
barrier gate shown on the left, which only depletes the upper 
well. The red arrow shows the direction of the externally ap¬ 
plied magnetic field B , which is in the plane of the wells and 
perpendicular to the wires. 


symmetrically in the AlGaAs barriers above and below 
the pair of wells. The doping is separated from the wells 
by spacer layers. The spacing between the centres of the 
two quantum wells is nominally d = 32 nm but we find a 
value of d = 35 nm fits the data better, and this can be 
explained by the fact that the centres of the wavefunc- 
tions will be slightly further apart owing to the opposite 
electric fields in each well. 

The 2DEG in the lower (dark blue) layer is used as 
a controllable injector or collector of electrons for the 
ID system4£ The current I tunneling in the ^-direction 
between the layers is proportional to the convolution of 
the ID and 2D spectral functions (a pair of peaks at 
k x = dz&p broadened in k y by the ID confinement, and 
a circle, respectively). An in-plane magnetic field B ap- 



Figure 7. Intensity plot of dG/dB at low energies around 
the Fermi point /cf- Spin (S) and charge (C) dispersions are 
indicated by dashed lines. The dotted lines indicate the para¬ 
bolae expected in the non-interacting model. The finger-gate 
voltage Vf — —0.70 V and the temperature T ~ 300 mK. 

plied in the ^/-direction, perpendicular to the wires (show 
with a red arrow in Fig. ED, produces a Lorentz force that 
changes the longitudinal momentum k x acquired while 
tunneling between layers by A k = eBd/h , where e is the 
electronic charge. Thus B shifts the spectral functions in 
k x relative to each other, and so probes the momentum. 
One spectral function can also be shifted relative to the 
other in energy by applying a voltage V between the lay¬ 
ers, in order to probe the ID and 2D dispersion relations 
at different energies. The conductance G = dl/dV has a 
peak when sharp features in the spectral functions have 
significant overlap. At V = 0 this occurs when A k is 
equal to the sum or difference of the Fermi wavenumbers 
kp and &2 of the ID and 2D systems respectively, so there 
are two peaks for B > 0, at B± = ±&f|- By sweep¬ 

ing B and V one can map out the dispersion relation of 
states in each layer. The range of magnetic fields that 
we apply to the system is still within the regime of Pauli 
paramagnetism for the electron densities in our samples. 

VII. LOW ENERGY 

First, we measure the tunnelling conductance G = 
dl/dB in a small range of voltages and magnetic fields 
around V = 0 and B = B+ = 3.15 T that corresponds 
to a region on the momentum-energy plane around the 
Fermi point (s = fi,k = &f), see Fig. 0 Below the 
Fermi energy we observe splitting of the single-particle 
line into two with different dispersions—spin (S) and 
charge (C) separatio n 11 ! 12 —giving two different slopes v s 
and v c (black dashed lines in Fig. [7]). We assume that v s 
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Figure 8. Intensity plots of dG/dV (left column, in /xS/mV) and dG/dB (right column, in /xS/T) from below — kp to above 
3/cf and from ~ —2/x to ~ 2/x, for various finger-gate voltages: —0.60 V (a, b), —0.65 V (c, d), —0.68 V (e, f) and —0.70 V 
(g-j). The solid black lines map out the dispersion of the lower (2D) layer. The green solid line marks a-modes, thick and thin 
dashed green lines, plb and hlb modes, respectively, and dashed blue, higher-A; modes (as in Fig. [3j. Dot-dashed yellow (blue) 
and dotted white (cyan) lines show second and third ID subbands (2D dispersion measured by those subbands), respectively 
(though the third is empty, electrons can tunnel into it from the 2D layer and hence there are sometimes signs of its effects for 
V > 0, especially near B — 0). Dotted magenta and blue lines are ‘parasitic’ 2D dispersions of the two layers. The voltage on 
the gate over this region Vp = 0 V except for e, f (Vp = 0.2 V) and i, j (Vp — 0.3 V), which shifts the parabolae to the right 
without changing the signal from the ID wires. The lines have all been adjusted to take account of the capacitive coupling 
between the layers. Spin (S) and charge (C) modes are indicated with black dashed lines. T ~ 300 mK. See Table ITT1 for the 
densities and the ratio v c /v s for each gate voltage. 


is the same as for non-interacting electrons and so take it 
to be the gradient of the parabola at V = 0. We estim¬ 
ate v c from the positions of steepest gradient and hence 
obtain v s ~ 1.2 x 10 5 ms -1 and v c ~ 2.3 x 10 5 ms -1 at 
the finger-gate voltage Vp = —0.70 V. 

Theoretically, the low-energy physics of the interacting 


ID electrons is described well by a spinful generalisation 
of the Luttinger-liquid models Its excitations are collect¬ 
ive hydro dynamic-like modes that are split into charge- 
only and spin-only excitations. For any finite strength of 
the interactions between fermions the two types of modes 
have linear dispersions with different slopes v c and v s . In 
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Vf (V) 

^ 2 d( 10 15 m -2 ) 

ni D (10 7 m -1 ) 

/cF(10 7 m x ) 

Vc/Vs 

-0.60 

1.67 

5.68 

8.9 

1.5 

-0.65 

1.65 

4.99 

7.8 

1.6 

-0.68 

1.52 

4.79 

7.5 

1.5 

-0.70 

1.48 

4.60 

7.2 

1.8 


Table II. Densities of the 2D layer (ti 2 d) and of the ID wires 
(tud), the ID Fermi wavevector kp (all to about ±1%), and 
the ratio of the charge and the spin velocities at low energies 
(to about ±5%), extracted from the gradients of the S and C 
lines in Fig. [8] for different finger-gate voltages Vp. 

the absence of interactions the difference between the two 
velocities vanishes in accordance with the free-electron 
model, in which the spin degree of freedom does not af¬ 
fect the spectrum but results only in the double degener¬ 
acy of the fermionic states. Thus the ratio of v c /v s serves 
as a good measure of the interaction strength. Since the 
Coulomb interaction between electrons is repulsive the 
charge branch always has a steeper slope v c > v s (see 
Ref. [ 4 ). Thus the ratio varies from 1 for free to 00 for 
infinitely repelling particles. In our experiment we meas¬ 
ure the tunneling of electrons and observe two peaks that 
we attribute to the charge and the spin dispersions. The 
pair of velocities above gives a large v c /v s « 1.8 =b 0.1 
(for Vp = —0.70 V), confirming that our system is in the 
strongly interacting regime. 


VIII. HIGH ENERGY 

Now we extend the ranges of the voltage and magnetic 
field measuring the tunneling conductance G across the 
double quantum well in Fig. [ 6 ] accessing a large portion 
of the ID spectral function from below — kp to 3kp and 
from —2/x to 2/x, see Fig. [ 8 j There is an unavoidable 
‘parasitic’ (‘p’) tunneling from narrow 2D regions (light 
blue strips in Fig. [ 6 ]) that connect the wires to the injector 
constriction. This superimposes a set of parabolic 2D-2D 
dispersions on top of the 1D-2D signal, which are marked 
by magenta and blue dotted lines in Fig. [ 8 ] 

Apart from the parasitic and the 2 D dispersion signals, 
we observe only a single ID parabola away from B — 0, 
marked by the solid green line in Fig. [ 8 j It extends from 
the spin-excitation branch at low energy and the position 
of its minimum multiplied by the electronic charge e gives 
the ID chemical potential /x ~ 4meV. The B _ and B+ 
crossings with the line V = 0, corresponding to momenta 
—kp and kp , give the ID Fermi momentum kp ~ 8 x 
10 7 m _1 . All other edges of the ID spectral function are 
constructed by mirroring and translation of the hole part 
of the observable ID dispersion, dashed green and blue 
lines in Fig. [ 8 j 

For positive voltages in the region just above the higher 
V = 0 crossing point (B + , which corresponds to kp) we 
observe a distinctive feature: the ID peak broadens, in- 
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Figure 9. Left column: intensity plots of dG/dV (in /xS/mV), 
for various finger-gate voltages and samples: (a, b) Vf = 
—0.68 V, Vp = 0.2 V, from Fig. [8^, sample A, which had 
10 /im-long wires (T ~ 300 mK); (c, d) Vp — —0.70 V, Vp = 
0.3 V, from Fig. [8j, sample A; (e, f) a similar single-subband 
result from sample B (18 /xm-long wires, T < 100 mK). The 
replica feature just above kp appears as a pale triangle (slowly 
varying G) between the two green curves, after a red region 
(sharp rise in G). The replica feature for sample B is some¬ 
what weaker than that for sample A, in line with the wire- 
length dependence predicted in this paper. Right column: G 
vs V at various fields B from 3 to 4.8 T for the data in the 
matching plots in the left column; ‘+’ and ‘x’ symbols on 
each curve indicate, respectively, the voltages corresponding 
to the dashed and solid (plb and pla{l)) green lines in the 
left column (and in Fig. [8|), showing the enhanced conduct¬ 
ance between the two. 


stead of just continuing along the non-interacting para¬ 
bola, with one boundary following the parabola (pla(l)) 
and the other bending around, analogous to the replica 
pib. This is visible in the conductance, but is most easily 
seen in the differentials, particularly dG/dV (left column 
of Fig. [ 8 ]). The broadening is observed at temperatures 
from 100 mK up to at least 1.5 K, and in samples with 
different wire designs (with or without air bridges) and 
lengths: in Fig. [9j dG/dV is shown in detail for the 
broadened ‘replica’ region for the 10 /xm wires already 
presented (a-d), and for another sample with wires 18 /xm 
long (e, f). G is plotted in Fig. (9b, d and f on cuts along 
the V axis of the corresponding plots in the left column 
at various fields B from B+ to 4.8 T—between the ‘+’ 
and ‘x’ symbols on each curve is a region of enhanced 
conductance characteristic of the replica pib. 

Filling of the second ID subband changes significantly 
the screening radius for the Coulomb interaction poten- 
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tial in the first ID subband. This is manifested by a 
change of the ratio v c /v 8 when the occupation of the 
second subband is changed by varying voltage of the fin¬ 
ger gates Vp in Figs. [E^ceg, see Table HT1 The ratio v c /v s 
is a measure of interaction energy. Thus, the finger gates 
give a degree of experimental control over the interactions 
within our design of the ID system. We use the max¬ 
imum change of the ratio v c /v s for different finger gate 
voltages to estimate the relative change of the interac¬ 
tion strength as (ma,x(v c /v s ) — min(v c /v s )) /min(v c /v s ) 
obtaining a change of about 20 %. 

It also has to be noted that the ‘replica’ is visible even 
when a second subband is present in the ID wires, see 
Fig. [ 8 ti-f. In a and b it appears to go 25-30% higher in 
voltage than expected for a precise copy of the usual ID 
parabola (even allowing for capacitive correction) due to 
a contribution of the second subband, which we do not 
analyse in detail here. 

At even higher magnetic fields the plb line passes a ‘p’ 
parabola. Figs. [9^i and c (and the corresponding cuts b 
and d) show the replica feature for two different positions 
of the ‘p’ parabolae using a gate above most of the ‘p’ 
region, showing that the replica feature is independent of 
the ‘p’ tunneling. The amplitude of the feature dies away 
rapidly, and beyond the ‘p’ parabolae, we have measured 
up to 8 T with high sensitivity, but find no sign of any 
feature that can be distinguished from the decaying tails 
of the other features. 

In the range of fields where the plb feature is observed 
its strength decreases as the B field increases away from 
the crossing point analogously to the power-law for spin¬ 
less fermions in Table HI On general grounds it is natural 
to expect that divergence of the spectral weight of a b- 
mode toward an a-mode is a general feature, but there 
is no known method for performing a microscopic calcu¬ 
lation in the spinful case. A similar feature should mark 
the hOb (r) mode (see Fig. [3] and Table H]) for negative 
voltages and for the magnetic field just below the cross¬ 
ing point kp , but it would be very difficult to resolve due 
to the overlaying spin and charge lines. 

Making an analogy with the microscopic theory for 
spinless fermions in the first part of this paper, we estim¬ 
ate the ratio of signals around different spectral edges 
using the ID Fermi wavelength, Ap ~ 130 nm for our 
samples, as the short-range scale 7 Z. The signal from the 
principal parabola, see Fig. [ 8 ) 3 , gives the amplitude of 
the a-mode as G a ~ 5 pS. Then the amplitude of the 
signal from the second (third)-level excitations is pre¬ 
dicted to be smaller by a factor of more than A p/L 2 ~ 
2 x 10 -4 (Ap/L 4 = 3 x 10 -8 ), where the length of a 
wire is L = 10 pm. These values G a A|/L 2 ^ 10 -3 /iS 
(G a Ap/L 4 ^ 10- 7 /iS) are at least two orders of mag¬ 
nitude smaller than the background and noise levels of 
our experiment G no i S e ^ 10 -2 /iS, which places an upper 
limit on the amplitude of any replica away from kp. Thus, 
our observations are consistent with the mode hierarchy 
picture for fermions. 

In an effort to quantify the decay of the replica fea- 
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Figure 10. (a,c) The conductance for Vf = —0.70 and 

—0.681/, respectively, after subtraction of an idealised land¬ 
scape made up of fits or estimates of the non-interacting 1D- 
2D and ‘p’ parabolae (see text). The plb replica is seen clearly 
as the red region of enhanced conductance, (b) The conduct¬ 
ance along the plb replica parabola, for the data in (a) (green 
crosses). The conductance on plb has a large contribution 
from the ‘p’ region (the line in (a) marked with blue dots, 
which is blurred to the left by multiple copies at slightly dif¬ 
ferent positions). In order to correct for this contribution, the 
conductance along a matching parabola shifted along the dot¬ 
ted ‘p’ line in (a) (shown as a dashed magenta line there), is 
subtracted from the plb data. This yields the points marked 
with blue circles, which appears to be non-zero because of the 
enhancement at plb. The amplitude decays rapidly. There are 
many uncertainties in the fitting of the other peaks, but the 
replica appears clearly and the decay of the conductance is 
consistent with an inverse-square power law G oc (k — kp)~ a 
(labelled a — 2 ), which is the behaviour predicted by the 
theory for k > kp + 7 where 7 <C kp (see Table [I]) . (d) The 
plb conductance enhancement as shown with circles in (c). 
Three different methods of fitting the background and the ID 
and 2D peaks are compared for each of two gate voltages as 
shown. The curves are offset vertically for clarity. The lines 
marked with values of a are guides to the eye. The data are 
all consistent with a = 2 d= 1 . 


ture we have fitted the gradual background fall in con¬ 
ductance and the non-interacting ID and 2D peaks (solid 
green and blue lines in Figs. EH2) with a Gaussian and/or 
Lorentzian functions of 5, at each value of V >0. The 
fitting parameters are then fitted to smooth functions in 
order to represent the general behaviour of the peaks as 
a function of V. This idealised landscape is then sub¬ 
tracted from the data, see Fig. flOh . and the ‘replica’ is 
then fairly easily observed in the remaining conductance. 
A copy of a nearby region along the ‘p’ curve is then 
subtracted too, as an approximation to the rather dif¬ 
fuse signal arising from the main ‘p’ peak and smaller 
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versions of it at slightly different densities. This also 
reduces errors in the peak and background fitting used 
in (a). We then plot the conductance along the ex¬ 
pected parabola (dashed line in (a)) as a function of 
(k - k F )/k F = (B - B+)/((B+ - B-)/ 2). This data 
is shown as circles in (c), where ah the other contribu¬ 
tions to the conductance along the same parabola are 
shown. Here, B+ = 3.17 T and k F = 7.2 x 10 7 m _1 . It is 
very hard to be sure that this procedure is reliable due 
to significant error bars imposed by contributions from 
the various other peaks, but it is clear that the replica 
feature dies away rapidly as a function of k — k F , and it 
is consistent with the 1 /(k — k F ) 1 2 3 4 5 6 law predicted for plb 
in Table I] for k — k F 7. Though the overall prefactor 
is unknown theoretically in the spinful case, this singu¬ 
lar power law may overcome the reduction factor 7 Z 2 / L 2 
close to k F . 

IX. CONCLUSIONS 

In this work, we have shown that a hierarchy of modes 
emerges in systems of interacting fermions in one di¬ 
mension at high energy controlled by the system length, 
in marked contrast to the well-known fermionic quasi¬ 
particles of a Fermi liquid and hydrodynamic modes of a 
Luttinger liquid at low energy. 

We have obtained theoretically the dynamic response 
functions for a model of spinless fermions with short- 
range interactions using the exact diagonalisation meth¬ 
ods of the Bethe ansatz for the spectrum and the form 
factors of the system. Analysing the spectral function in 
detail, we have found that the first-level (strongest) mode 
in long systems has a parabolic dispersion, like that of a 
renormalised free particle. The second-level excitations 
produce a singular power-law line shape for the first-level 
mode and different kinds of power-law behaviour at the 
spectral edges. Evaluating the form factor necessary for 
the dynamical structure factor we have shown that it has 
the same general form as the form factor of the spectral 
function, manifesting the same hierarchy of modes. 

Using the same many-body matrix elements obtained 
microscopically, we have also calculated the local density 
of states. It provides a more convenient way to ana¬ 


lyse how the hierarchy at high energy changes into the 
hydrodynamic modes of the Luttinger liquid at low ener¬ 
gies. We have shown, via a full Bethe-ansatz calculation, 
that the LDOS is suppressed at the Fermi energy in a 
power-law fashion in full accord with the prediction of 
the Tomonaga-Luttinger model. Away from the Fermi 
point, where the Lorentz invariance of the linear disper¬ 
sion is reduced to Galilean by the parabolicity of the 
spectrum, the LDOS is dominated by the first (leading) 
level of the hierarchy. We have demonstrated that the 
transition from one regime to another is a smooth cross¬ 
over. 

We measure momentum-resolved tunnelling con¬ 
ductance in one-dimensional wires formed in the 
GaAs/AlGaAs double-well heterostructure by an array 
of finger gates. In this set-up we probe the spectral func¬ 
tion of unpolarised electrons (spinful fermions) and find 
a pronounced spin-charge separation at low energy with 
a ratio of the spin and the charge velocities up to 1.8, 
which confirms that our system is in the strongly inter¬ 
acting regime. By varying the gate voltage that controls 
the width of our ID wires, we demonstrate control of 
the interaction strength of about 20%; the deeper confin¬ 
ing potential of the wires populates higher ID subbands 
as well which in turn screens stronger Coulomb interac¬ 
tions in the principal ID band reducing the interaction 
strength. In 10 /im-long wires we find a clear feature re¬ 
sembling the second-level excitations, which dies away 
rapidly at high momentum. A qualitative fit shows that 
the feature decays in a fashion that is consistent with the 
power-law prediction in this paper for spinless electrons. 

Thus we have shown that the hierarchy is apparently 
a generic phenomenon at least for one- and two-point 
correlation functions of fermions without spin, and for a 
transport experiment for fermions with spin. 
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Appendix A: Eigenvalue equation in the algebraic framework 


The eigenvalue of the transfer matrix r(u) in Eq. ([T9l) can be evaluated using the commutation relations in Eqs. 
(116I17D . Let T be a Bethe state in the algebraic representation of Eq. (|2T1) . The results of acting with A (u) and 
D ( u ) operators on the state 4/ are obtained by commuting these operators from left to right though the product of 
C ( Uj ) in Eq. ([2T1) and then by using the their vacuum eigenvalues in Eq. (l20l) . 

Let us consider the case of N = 2 and the operator A ( u ) first. Commuting once by means of Eq. m gives 


A(«)C(«2)C(ui)|0) = 
Applying Eq. (11711) the second time gives 
A(u)C(u 2 )C( Ul ) |0)= ( 


hM) c MA{u) - <“> A w ' c <“>> i°> • < A1 > 

C ( U\ — u) 1 


+ 


b {u 2 — u)b (u\ — u) 

1 C{U2~ U) 

b (ui - U 2 ) b {u 2 - u) 


C (u 2 ) C (ui) a (u) 


b (ui — u) b (u 2 — u ) 


C (u 2 ) C (u) a (ui) |0> (A2) 


C (■u) C(u x )a M + -, C («) C (u 2 ) a ( Ul ) ) |0> 

b [ui - U2) b [U2 - u) ' 


where the vacuum eigenvalue of A (u), A (u) |0) = a (u) |0), was substituted explicitly. The second terms in the first 
and the second lines of the above equation have the same operator form but different scalar factors. Summation of 
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the two scalar factor, using the explicit form of b(u) and c(u) from Eq. ([T2l) , yields 


a (ui) 


c(ui — U2) , v c(ui — u) 

C[U 2 - U) - 


a (u\) c (ui — u) 


b(u2~u) [b(ui—U2) b{ui—u)\ b (U2 — ui) b (ui — u) 

Thus the four terms can be rewritten as only three terms, 

A (") c M c («,> io) =(«m n - g ° (“j > | c ( “> ,_n a _ tj > ) |o> - 


i=l#J 


Extension of the same procedure for TV > 2 gives 

at f n 1 

a(«) n c ( u j) i°) = °( y ) n uz ~a c 

3 = 1 \ i=i 1 J J 


AT 




AT 


3=1 ' Z=l#j 

Commuting of the operator D (u) is done in the same way using Eq. (ED and yields 


b(ui - Uj ) 


(A3) 


(A4) 


<?(««) 10) (A5) 


jv / JV _ JV C _'\ JV _1 \ 

CW nc ( ,,)| 0 )= + JI 10 ) 


(A6) 


Thus a Bethe state in Eq. (ED parametrised by an arbitrary set of Uj is not an eigenstate of the transfer matrix 
since acting of the operator r(u) on such a state does not only result in multiplying by a scalar but also generates 
many different states: the second terms in Eq. (1A5D and in Eq. (1A6D . However, the coefficients in front of each of 
these extra states can be made zero by a specific choice of Uj , 


N 

a{uj) H 


1 

b (ui - Uj) 


N 


d ( Uj ) n 

*=i#j 


-1 

b(uj - ui) 


= 0 . 


(A7) 


When a set of uj satisfies the system of equations above, the corresponding Bethe state is an eigenstate of the transfer 
matrix, t(u)\ P = 'T(u)'f? , with the eigenvalue T given by the first terms in Eq. (IA5II and in Eq. (IA6I) . 

N I N _i 

T(u) = a(«)U^^-d(u)j\^^ y (AS) 

Eq. (|A7D is the set of Bethe equations in the algebraic representation, Eq. (|22l) of the main part of the text, and Eq. 
(|A8I) gives the eigenvalue of the transfer matrix, Eq. ([24l) of the main part of the text. 


Appendix B: Calculation of averages of the local density operator p (0) 


The calculation of the average of the local density operator p (0) is done in the same way as for the field operator 
(0) in Subsection IVc. We start from the lattice model in Eq. 0 and the corresponding construction of the 
algebraic Bethe ansatz in Subsection IVa. 

The local density operator can be represented in terms of A and D operators a a 17 i 18 


pi = -D 



(Bl) 


The action of the second factor in the above expression on an eigenstate |u) just gives a phase factor - see an 
explanation after Eq. m - that we will ignore since we are interested in the modulus squared of this form factor. 
Then commuting the operator D of the first factor in the equation above through all C operators of the eigenstate 
|u) - in the form in Eq. (ED - gives the result in Eq. (|A6D . 
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The scalar product of Eq. (ESI, where the auxiliary parameter u is set to itt /2 — 77, with another eigenstate (v| 
gives 


N 


<v|pi|u) = (~1) N J 


N -TT COSh (Uj - 7 ?) 

cosh (uj + 77) 
1=1 v J ' 


(u|v) 


tv sinh 2 rj sinh ( u 5 — ui + 2 rj) 


1 if MllllZ// TT 

“ cosh (7/6 + 77) & sinh (7/5 - ui) 


ITT 


Ub- 1, y - 77,1/5+ilv ) , (B 2 ) 


where (7/5-1, ^ — 77,7/5+1! is a Bethe state which is constructed from the eigenstate u by replacing b th quasimomenta 
with in /2 — 77. Note that the properties (v|u) = (u|v), where Uj satisfy the Bethe equations and Vj is an arbitrary 


set of quasimomenta j 17 i 18 was used. 


The scalar product in the first line in Eq. ( 1 B 2 I is given by Eqs. (EH HEI) where u and v are swapped. Substitution 
of u = 7/5-1, — 77,7/5+1 in the same expressions for the scalar products in the second line of Eq. ( 1 B 2 D yields 


77 T 

Ub-i, y - V,u b+ i\v 


sinh^ ( 2 rj) cosh (vj — rj) det T^ 


II j<i sinh ( vj - Vi) sinh ( Uj - u») (-l) b 1 fl^Li cosh ( u j + V) 


(B 3 ) 


where all matrix elements of T (l>l are 
N 


T (b) , = TT 

ab ' sinh i u b' -Ut- 2 r?) sinh (u b > - v a ) * 


sinh (7/5/ — ui + 2 rj) 


1 


N 


sinh (7/5/ — Vj — 2 rf) 


N 


sinh (7/5/ - v a ) 


sinh (7/5/ — Vj + 2 rf) (B 4 ) 


j= l#a 


for b' ^ b 


/fiiP) _ 

± ab' 


1 

cosh (v a + 77) cosh (v a - 77) 


(B 5 ) 


for b' = b. 

After pulling a common factor out of the brackets in Eq. (jB 2 D and absorbing the 6 -dependent prefactors in front 
of the determinants in the second line of Eq. (|B 2 D into the b th columns of the matrices under the determinants, the 
form factor in Eq. (|B 2 D reads as 


( v l/h|u) 


N 


n 


cosh (Uj — 77) sinh^ ( 2 rj) 

cosh (Uj + 77) Yl j<:i sinh {vj - Vi) jKi sinh (uj - Ui) 


det T 


n . ' 

■^2 det f (6) 
5=1 


(B 6 ) 


Here the matrix elements of T^ b \ which are obtain by multiplying by the corresponding scalars, are 


N 

«= n 


sinh (7/5/ — ui + 2 rf) 


N 


11 n / v 7 -r TT sinh (7/5/ — Vj — 277 ) 

j sinh (7/5/ - ui- 277 ) sinh (7/5/ - v a ) 


N 


sinh (u b > - v a ) \ 

3=1 /a 


sinh (7/5/ — Vj + 277 ) (B7) 


for b' 7^ b 


rp(P) 

1 ab' 
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(-1)^ J4 sinh (uv - ui + 2r?) JJ 

i=1 


cosh (vj — 77) sinh (277) 

cosh (Uj — 77) cosh (u a + 77) cosh (i; a — 77) 


(B 8 ) 


for b' = b. Note that T^y = Tab' for b 7^ b. 















23 


Finally, the summation over b in Eq. (1B6D can be evaluated using a general matrix identity: det T + Ylb =1 det = 
det (i 1 + B^j where T^) is obtained from the matrix T by substituting b th column from matrix B. After constructing 

the matrix B out matrix elements from Eq. ()B 8 D and performing the summation over b in Eq. (|B 6 |) . the form 
factor reads 


N 


<v|pi|u> = J 


cosh (uj — rf) 


sinh^ (2rf) 


L = i c ° sh («i+ v) n j<i sinh & - n i<4 sinh ( u j - u i ) 


det if 


(B9) 


where the matrix elements of K are 
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K, 


ab 
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sinh (Ub — ui + 277 ) 
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N 


l=1 ^ b sinh ( u b — ui- 2 rj) sinh (u b - v a ) j=1 _^ a 


sinh (u b — Vj — 2rf) — 


sinh (ub - v a ) 


sinh {u b — Vj + 2rf) 


j = l/a 


N 


N 


(- 1 )^ tl sinh (u b - ui + 2rj) J 


cosh (vj — rf) 


sinh ( 277 ) 




^ cosh (uj — rf) cosh (v a + rf) cosh (v a — rf) 


(BIO) 


Now we evaluate the long wavelength limit for the result above. Applying the inversion mapping from the algebraic 
to the coordinate representation from Eq. to the matrix elements in Eq. (jBlOD and expanding the result up to 
the leading order fcj, kj <C 1 we obtain, 


K ab = (-if- 1 2mU (( mU) 2 - l)' 


E N h v — h u — h v -L h u 

j =1 2 ^j= 1 ^ 


+ 2 (mU + 1) fmUf - l) ' 


(Bll) 


Repeating the same procedure for the prefactor in Eq. (|B9D and pulling a common scalar factor out of the matrix 
elements under the determinant we obtain 


(v|p(0)|u) = (-l) Jv2 2 w 


(mU) 


N 


(mU — 1) N det A: 


(mU + 1 ) 


N 


(B12) 


where the matrix elements of /C are 


and A P = E^i K ~ T,f=i K- 


Kab = k l ~ k l ~ AP + mU + 1 


ibti _ ib'U 


mU 


(B13) 


j = l ,XJ j = l 3 • 

Evaluating the determinant we obtain 


det /C = 


mu +1 (a pf n,<, (k? - fcp n,<, {k - *?) 


mU 


n i,j {k - k?) 


(B14) 


This formula can be proved by induction analogously to the proof of Eq. (|43|) . 

The form factor appearing in the dynamical structure factor in Eq. m is a modulus squared of Eq. (IBT 211 . 

Normalising the initial and the finial states using Eq. (1381) as (f\p (0) |0 )| 2 = |(k/|/?(0) k°) | 2 1 (k°|k 0 ) 1 

we obtain 


/ tt \2N—2 P^ N 

K/Ip(o)|o)| 2 = (mU) 




(mC/+1)2JV_2 


(B15) 


where Pq = 0 for the ground state. Eq. (jB15[) is Eq. (1471) in the main part of the text with Z = 
mU/ (mU + 1 )/(£- NmU/ (1 + mU)). 

Note that when the final state becomes the ground state, k/ = k°, Eq. ()B15D is divergent. In this case the matrix 
element is evaluated using only translational symmetry and the definition of the number of particles operator as 


AT 2 

l(/|p(o)|o )| 2 = -^, 


(B16) 


which also follows directly from Eqs. (|B9i IB 10 D , by taking the limit v —>> u. 





















